A study of the association between single nucleotide polymorphisms of the endoplasmic reticulum aminopeptidase 2 (ERAP2) gene and the risk of ankylosing spondylitis in Egyptians

Background Ankylosing spondylitis (AS) is often regarded as the prototypical manifestation of spondylo-arthropathies that prevalently involves the axial skeleton with the potential attribution of ERAP2 polymorphisms to AS predisposition. The purpose of this study was to determine the genetic association between ERAP2 gene rs2910686, and rs2248374 single nucleotide polymorphisms (SNPs) and the risk of ankylosing spondylitis in the Egyptian population. Methods and results A cross-sectional work involved 200 individuals: 100 AS individuals diagnosed based on modified New York criteria in 1984 with 100 healthy controls matched in age and gender. The study included a comprehensive evaluation of historical data, clinical examinations, and evaluation of the activity of the disease using the Bath Ankylosing Spondylitis Disease Activity Index (BASDAI). A comprehensive laboratory and radiological evaluation were conducted, accompanied by an assessment and genotyping of the ERAP2 gene variants rs2248374 and rs2910686. This genotyping was performed utilizing a real-time allelic discrimination methodology.Highly statistically substantial variations existed among the AS patients and the healthy control group regarding rs2910686 and rs2248374 alleles. There was a statistically significant difference between rs2910686 and rs2248374 regarding BASDAI, BASFI, mSASSS, ASQoL, V.A.S, E.S.R, and BASMI in the active AS group. Conclusions ERAP2 gene SNPs have been identified as valuable diagnostic biomarkers for AS patients in the Egyptian population being a sensitive and non-invasive approach for AS diagnosis especially rs2910686. Highly statistically significant variations existed among the AS patients and the healthy control group regarding rs2910686 alleles and genotypes.Further research is recommended to explore the potential therapeutic implications of these SNPs.


Introduction
Ankylosing spondylitis (AS) is widely recognized as the prototypical variant of spondyloarthropathies.It is a persistent autoimmune-inflammatory condition characterized by progressive arthritis that primarily affects young adults.The disease predominantly involves the axial joints including the spine and sacroiliac joints, and is often accompanied by additional manifestations outside the joints, such as psoriasis, uveitis, and inflammatory bowel diseases [1] .
Because of the lack of precise understanding of the pathophysiology, it is hypothesized that the clinical progression of AS starts with an initial phase of inflammation, subsequently leading to the development of the formation of new bones that triggers localized osteitis, erosion of cartilage, destruction of bone, and eventual ankylosis.Sacroiliitis and syndesmophytes, known as radiographic AS, may be identified with the use of conventional radiography.However, the early diagnosis of these conditions can be facilitated by using MRI [2].
The human leukocyte antigen (HLA-B27) remains the most effective biomarker for diagnosing AS, whereas C-reactive protein (CRP) serves as the optimal marker for evaluating the extent of the disease, assessing therapy effectiveness, and monitoring structural development.However, it should be noted that HLAB-27 accounts for just 30% of the genetic variables associated with AS, suggesting the presence of other genetic abnormalities implicated in the development of AS [3].
The enzyme known as endoplasmic reticulum aminopeptidase 2 (ERAP2) is situated inside the endoplasmic reticulum.It is classified as a member of the zinc-containing metallopeptidase family, and its corresponding gene is positioned on chromosome 5q15.The involvement of ERAP2 in the process of peptide trimming via the major histocompatibility complex class I throughout the antigen presentation pathway has been observed.In contrast to ERAP1, there is a limitation of evidence on the association between ERAP2 polymorphisms and AS susceptibility [4].Several genetic variants within the ERAP2 gene have been associated with modifications in the structure and functionality of the protein.The SNP rs2248374 in the ERAP2 gene has been shown to have a protective effect against AS.This SNP specifically affects the splicing location in exon 10 of the gene, resulting in the production of a longer exon 10 transcript [5].
ERAP2 instead has evolved under balancing selection that maintains two haplotypes, one of which undergoes RNA splicing leading to nonsense-mediated decay and loss of protein.Hence, likewise in rodents, wherein the ERAP2 gene is missing, about a quarter of the human population does not express ERAP2 [6].
The rs2248374 SNP has been previously associated with a reduced risk of AS.This SNP is a functionally significant variation located in the ERAP2 gene, that plays a crucial role in modifying the functional velocity and specificity of ERAP2 functioning during the process of trimming peptides.The rs2910686 SNP located inside the ERAP2 gene has been identified as a gene of functional significance, and its association with AS susceptibility has been established.The association between SNPs in the ERAP1 gene and susceptibility to AS has been reported [7].
The ERAP2 gene variant rs2549782 SNP exhibits linkage disequilibrium (LD) with many other ERAP2 SNPs, such as rs2548538, rs2287988, rs1056893, and rs2248374.These SNPs serve as marker variants that form haplotypes A and B, which are correlated with the protein expression of ERAP2.Furthermore, the ERAP2 gene variant rs17408150 results in a nucleotide substitution from T to A at codon 669, leading to the amino acid change from leucine to glutamine (p.Leu-669Gln).This alteration has been seen to have a significant impact on the functionality of the ERAP2 enzyme [5].
The objective of this work is to investigate the potential genetic correlation between two SNPs, namely rs2910686 and rs2248374, in the ERAP2 gene and the susceptibility to AS within Egyptian patients with AS.To the best of our knowledge, this is the first attempt to explore this link in this specific group.Additionally, an investigation was conducted to explore the potential involvement of these SNPs in regulating inflammatory markers.

Patients and methods
Based on a previous recent study (Ebrazeh et al., 2021), found that in the instance of rs17408150, the C allele had a substantial association with a greater risk of AS in the HLA-B27-positive individuals (OR1.39,95%CI1.06-1.81,P = 0.013).The last sample size calculated utilizing statistics and sample size version 6 is 200 subjects, divided into two equal groups.The power of the work is 80% and the confidence level is 95% [8].
This cross-sectional work involved 200 individuals: 100 AS individuals diagnosed based on modified New York criteria in 1984 [9] recruited from the rheumatology and clinical immunology outpatient clinic, the outpatient clinic of the Rheumatology, Physical Medicine, and Rehabilitation Department, in collaboration with the Biochemistry and Molecular Biology Department, in the duration between January 2022 and July 2022 with 100 healthy controls matched in sex and age.
The present work received approval from the institutional research ethics council at our university, with the assigned IRB number 1/2022INTPH6-1.A comprehensive informed consent was collected from all participants included in the present study.
Individuals with various autoimmune illnesses, active infection, lymphoproliferative disorders, cancer, myocardial infarction, heart failure, ischemic injury, and thrombus formation, as well as patients with diabetic nephropathy, were not included in this study.
Demographic information was collected for all participants.The length of disease, specific behaviors, including smoking, clinical evaluation, and past therapeutic interventions were obtained from individuals diagnosed with AS.The assessment of the quality of life was conducted utilizing the AS quality-of-life questionnaire (ASQoL) [10], in addition to pain measurement recorded on a 10 cm visual analog scale (VAS) [11].
The evaluation of the disease's activity was conducted by using the BASDAI, in conjunction with measurements of CRP levels in (mg/l) and ESR in mm/h [12].The assessment of mobility and functional constraints was conducted using the Bath AS Metrology Index (BASMI) [13] and the Bath Ankylosing Spondylitis Functional Index (BASFI), respectively [14].
The participants were categorized into two distinct groups using the criteria of ESR, CRP, and BASDAI.The first group consisted of active AS individuals who exhibited a CRP level of more than 8 mg/l and/or a BASDAI score of 4 or higher, together with an ESR above 20 mm/h while the remainder of the patients were categorized as inactive AS instances.
Comprehensive laboratory analyses were conducted, which included genetic variables such as HLA-B27, antigens, inflammatory indicators (CRP and ESR) measured by the Westergren technique, and 25-hydroxyvitamin D [25(OH) D] concentrations assessed via the enzyme immunoassay method utilizing a homogeneous enzyme linked to vitamin-D binding protein.
Lateral radiographs were conducted on the cervical and lumbar spines, and the evaluation of structural harm was performed using the modified Stoke Ankylosing Spondylitis Spinal Score (mSASSS).This scoring system involves assessing the existence or lack of syndesmophytes, along the anterior margins of the lower border of C2 to the upper border of Th1, as well as from the lower border of Th12 to the upper border of S1.Each syndesmophyte is assigned a grade ranging from 0 to 3 points, with 0 indicating normal findings, 1 indicating erosion, sclerosis, or squaring, 2 indicating the existence of a syndesmophyte, and 3 indicating the formation of a bone bridge.The total score ranges from 0 to 72 [15].
Blood specimens: 6 ml of blood samples were taken under strict aseptic conditions from each person using clean venipuncture.Each sample was divided into three tubes as follows: Tube 1: EDTA-contained sterile tube kept at -20 C for further DNA extraction for HLAB27 analysis and genotyping of rs2248374 and rs2910686 of ERAP2 gene, Tube 2: Sodium citrate sterile tube presents at room temperature for measurement of ESR, Tube 3: sterile plain tube that allowed to clot at 37° C, and used for measurement of CRP and vitamin D after serum separation by centrifugation.
DNA was extracted from whole blood by a GeneJET Genomic DNA Purification Kit (Thermo Scientific, Lithuania.cat# K0721) following the manufacturer's protocol.DNA concentration, quality, and purity were assessed using a Nanophotometer N-60 (Implen, Germany).
HLAB27 analysis was done using the HLA-B27 REAL-TIME PCR Genotyping Kit with a PCR mix that includes the internal control (IC B27).ICB 27 serves as sample intake control and allows to evaluation of the quantity of genomic DNA.It is needed for assurance of PCR quality and sufficiency of input DNA.The interpretation of the result for each sample is carried out automatically concerning the Threshold cycle (CT) for Fam (specific dye label) and Hex (internal control dye label) channels.
ERAP2 gene rs2248374 and rs2910686 SNPs genotyping was done using an allelic discrimination assay by realtime PCR using a TaqMan probe (Applied Biosystems, Foster City, USA) with TaqMan SNP Genotyping Assays.Master Mix II (2x), primers, and probes were supplied by Thermo Fisher Scientific and Catalog number (Catalog No D3024).(Fig. 1a and b, respectively)

Detection of ERAP2 gene rs2248374 and rs2910686 SNPs
Thermal cycling conditions: Following the PCR amplification process, the endpoint plate was subjected to analysis utilizing the Applied Biosystems Real-Time PCR System.The Sequence Detection System Software utilizes the fluorescence measures obtained throughout the plate read to provide a graphical representation of fluorescence (Rn) values, which are derived from the signals detected in each well.The fluorescence signals that have been displayed provide information on the presence of certain alleles in all samples.

Statistical analysis of the data
The data were inputted into the computer and then analyzed utilizing the IBM SPSS software program, version 20.0.(IBM Corp, Armonk, NY).The representation of categorical data included the use of numerical values and percentages.The chi-square test was used to conduct a comparative analysis among two distinct groups.In cases where the number of controls, matched in terms of age and sex, also participated in the study.

Demographic and laboratory data of patients and healthy controls
There were significant differences among active ankylosing spondylitis (group I) and healthy control (group III) regarding smoking habit (P < 0.001).additionally, a significant increase was found in active ankylosing spondylitis, inactive ankylosing spondylitis (group II), and healthy control regarding ESR, CRP, and vitamin D levels, but there were non-significant differences between active ankylosing spondylitis, inactive ankylosing spondylitis and healthy control regarding age, sex, occupation, disease duration, and drug taking (Table 1).

Univariate regression analysis for alleles of rs2248374 and rs2910686 between cases and controls
rs2248374 SNP G allele was more frequent in the AS patient group (54%) compared to the healthy control group (49.5%), but the difference was not statistically significan.The frequency of the rs2910686 SNP C allele was significantly higher in the AS patient group compared to the healthy control group ( 94(47%) versus 50(25%), respectively) with OR predicted count of < 5 was seen in over 20% of the cells, the Fisher Exact or Monte Carlo adjustment test was used as an alternative method.The normality of the continuous data was assessed using the Kolmogorov-Smirnov test.The quantitative data were presented in terms of the range (minimum and maximum values), mean, standard deviation, and median for parameters that followed a normal distribution.The F-test, (ANOVA), has been employed to assess the differences across many groups.Additionally, a Post Hoc test, specifically the Tukey test, was conducted to do pairwise comparisons among the groups.In contrast, when dealing with quantitative parameters that do not follow a normal distribution, the Whitney test was employed for contrasting two groups, whilst the Kruskal-Wallis test was utilized to contrast more than two groups under study.Additionally, pairwise comparisons were conducted using the Post Hoc method, specifically Dunn's multiple comparisons test.The significance of the obtained results was judged at the 5% level.

Results
A total of one hundred individuals diagnosed with AS were categorized into two groups, namely active and inactive, using criteria such as ESR, CRP, and BASDAI.A similar mSASSS(P < 0.001), ASQoL(P < 0.001), V.A.S(P < 0.001) in individuals with active ankylosing spondylitis when contrasted to those with inactive ankylosing spondylitis.Also, there is an increased percentage of patients with clinical symptoms in the form of Knee arthritis, Enthesitis, and Uveitis with an x-ray indicating Ankylosed cervical spine in individuals with active ankylosing spondylitis when contrasted to inactive ankylosing spondylitis (Table 3).
Significant differences were also observed for the frequencies of the genotypes for this SNP between AS cases and controls.The number of carriers of the minor allele C (genotypes TC and CC) was significantly higher in patients compared to controls (72% versus 42%, P < 0.001) (Table 2).

Discussion
Ankylosing spondylitis (AS) is a chronic inflammatory autoimmune condition that is distinguished by the presence of pathologic new bone growth [16].The origin, pathophysiology, and diagnosis of AS continue to provide significant challenges in the field.The HLA-B27 antigens have been widely seen in individuals with AS, making it a significant diagnostic biomarker.However, it is important to note that the contribution of this factor to the total genetic susceptibility of AS is below 30%.In addition, it is important to note that X-rays are not capable of detecting first radiological alterations, while the use of MRI as a screening tool is hindered by its high cost and impractical handling.Hence, there is a critical requirement for the identification of new biomarkers to facilitate the screening and diagnosis of AS.
The accurate diagnostic and prognostic considerations of AS remain unaddressed medical needs.The prioritization of

Relation between (rs2248374) gene polymorphism and demographic, laboratory data, and clinical data of patients with inactive ankylosing spondylitis
There were no statistically significant differences among different parameters and the three types of genotypes (Table 6).The cohort of participants with active AS had a lower mean age (35.7 ± 7.9 years) compared to those with inactive AS.Furthermore, the active AS group had a longer average duration of illness (25 ± 2.4 years) and mostly consisted of male individuals (69%).This gender distribution aligns with the typical epidemiological trend of AS, whereby men are afflicted at a frequency of 2 to 4 times higher than females.Our findings are consistent with evidence from the literature that has been published and reflects the disease's youthful character, which makes it more aggressive in young people and has a peak age of onset between the 2nd and 3rd decades [15].Numerous investigations have implicated the male gender and extended illness duration as unfavorable prognostic indicators for ankylosing spondylitis [18].But care for individuals with AS is crucial due to the detrimental effects of this chronic inflammatory condition on young people, posing risks to their daily activities, productivity, and overall quality of life.This, in turn, leads to an enormous financial burden on a national scale [16,17].
The primary objective of this research is to investigate the potential genetic correlation between two specific SNPs, namely rs2910686 and rs2248374, inside the ERAP2 gene, and the susceptibility to AS.To the best of our knowledge, this study represents the first attempt to explore this link in a cohort of 100 individuals with AS within our nation.Additionally, an investigation was conducted to explore the potential involvement of these SNPs in regulating inflammatory markers.vitamin D may play a major role in the development and progression of the illness.This statement aligns with several recent studies that have provided evidence for the involvement of vitamin D in AS via influencing both the adaptive and innate immune systems, rather than only affecting calcium metabolism.These findings suggest the potential also it was shown that both males and females had functional incapacities and structural damage [18].
In the present study, it was observed that there was a deficiency of vitamin D in the active and inactive groups of individuals with AS, correspondingly.Notably, there was a substantial decrease in the number of active AS instances (16.05 ± 2.73), which suggests that inadequate levels of refuted the notion of vitamin D having any immunomodulatory impact on AS disease activity.HLA-B27 has traditionally served as the primary diagnostic marker for AS, despite the recognition that this class 1 human leukocyte antigen accounts for just 30% of the overall hereditary risk associated with AS [21].The research benefits of supplemental vitamin D for individuals with this particular condition [19,20].
There are conflicting literature results about the correlation between vitamin D and ankylosing spondylitis (AS), as shown by the research conducted by Zhao et al., [20], which  [19].The findings may validate the pursuit of novel diagnostic markers for AS that provide improved accuracy.
In the present work, we revealed that there were highly statistically significant differences among the active AS patients and inactive AS patients regarding BASDAI, BASFI, BASMI, mSASSS, ASQoL, and V.A.S. conducted by Neves et al., revealed that HLAB-27 had a high prevalence rate of 96.8% among a cohort consisting of individuals of Caucasian descent.The prevalence of HLAB-27 in AS sufferers was found to be 66% in this research.The observed disparity might perhaps be ascribed to the racial genetic variability within the groups under investigation    In their study, Londono et al. observed statistically significant variations among the patients under investigation in terms of their scores on BASDAI and BASFI, as well as the overall assessment of the disease's activity conducted by the physician [23].
This research provides evidence of substantial and statistically significant variations between the groups of Our findings were in line with Maghraoui et al., who found that statistical disparities were seen in symptomatic severity measures, including the Schöber test, BASMI, chest expansion, BASFI, BASDI, and BASG, as well as structural severity metrics, such as lumbar syndesmophytes score and BASRI when comparing patients with active AS to patients with inactive AS [22].In contrast to our findings, Su et al. found no association between SNPs of ERAP1/ERAP2/RUNX3 and susceptibility of AS with or without AAU.A case-control study between patients with human leucocyte antigen HLA-B27-positive and healthy controls also failed to demonstrate an association of the tested SNP with AS with or without AAU.Moreover, a meta-analysis showed that there was no association between rs30187, rs27037, rs27980, rs27434, and rs27582 in ERAP1 with AS in Chinese Han.Taken together, 17 SNPs in ERAP1/ERAP2 and RUNX3 genes did not confer disease susceptibility to AS in Chinese Han [27].
Gao et al. also did not reveal any evidence of the associations of rs2248374 and rs2549782 polymorphisms in the ERAP2 gene and susceptibility to AS.This is the only published meta-analysis of ERAP2 SNPs in AS including A total of six studies including 2774 AS patients and 4119 disease-free controls were eligible for this meta-analysis.Five studies reported rs2248374 polymorphism and three studies reported rs2549782 polymorphism.There were no statistically significant association with AS susceptibility: rs2248374, A vs. G, OR = 0.94, 95%CI 0.86-1.02,P = 0.14; rs2549782, T vs. G, OR = 1.03, 95%CI 0.95-1.12,P = 0.45 [28].
This research provides evidence of a statistically significant association between the rs2248374 genetic variant and variables related to occupation, drugs, smoking, and nonsteroidal anti-inflammatory drugs (NSAIDs).The findings of our study were corroborated by the research conducted by Bugaj et al., [29], which demonstrated a statistically significant association among rs2248374 and NSAIDs (p < 0.05).Furthermore, a strong correlation between rs2248374 and smoking was observed by Robinson et al., (p < 0.05) [24].
This research elucidated that a statistically significant disparity existed about rs2910686 concerning smoking, BAS-DAI, BASFI, and BASMI.The findings of our study were corroborated by Bugaj et al., [29], as they also observed significant relationships between the rs2910686 genetic variant and both the VAS and BASDAI scores.
This research provides evidence of a statistically significant disparity in the variables of rs2910686 concerning mSASSS, C.R.P, E.S.R, and Vitamin D.
In a similar vein, Bugaj et al., [29] discovered a noteworthy association between rs2910686, ESR, and CRP, with a statistical significance seen at a p-value of less than 0.05.Moreover, the study conducted by Robinson et al., [24] revealed a statistically significant correlation between rs2910686 and vitamin D levels (p < 0.05).
by univariate logistic regression, the rs2248374 SNP G allele and the rs2910686 SNP C allele were significantly upregulated in the AS patient group compared to the healthy control group with OR (LL -UL 95% C.I) of 1.198 (0.809-1.774) and 2.660 (1.741-4.064),respectively.Our findings individuals with AS and the control group in E.S.R and C.R.P.In a similar vein, Ebrazeh et al., [8] observed a statistically significant elevation (P < 0.0001) in the C-reactive protein level among those in the case group (3.14 ± 2.27) compared to those in the control group (1.27 ± 0.95).Furthermore, Robinson et al., [24] discovered a statistically significant elevation in levels of E.S.R and C.R.P in individuals with AS when contrasted with the control group (p < 0.001).
In a study conducted by Robinson et al., [24] in 2015, the authors examined specific subgroups of HLA-B27 positive AS-associated ERAP1 haplotypes, which were identified by the genetic markers rs30187 and rs10050860.The study also included HLA-B27 positive and unselected controls of European descent from the IGAS study.The researchers discovered a correlation between the ERAP2 polymorphism rs2248374 and AS in individuals who tested positive for HLA-B27.These investigations revealed a significant correlation between ERAP2 and HLA-B27-positive AS individuals.
The research conducted by Cortes et al., [25] in 2013 included a total of 10,619 patients diagnosed with AS and 15,145 control subjects of East Asian, European, and Latin American descent.The researchers identified a correlation between AS and two SNPs located in the ERAP2 gene, namely rs2549782 and rs2248374.This link was seen among individuals who tested positive for the HLA-B27 antigen.Upon doing a controlled analysis of the ERAP1 interaction with the ERAP2 polymorphism, a significant connection was identified between ERAP2 SNPs and AS instances.Specifically, the rs2910686 SNP exhibited an OR of 1.19 and a p-value of 0.213 × 10 − 2. The study conducted by Ebrazeh et al., [8] revealed a substantial statistical connection between the rs2910686 SNPs and the risk of AS.
In the present study, there were significant associations between ERAP2 polymorphism rs2248374, rs2910686, and risk of AS.A statistically significant difference was observed in the frequency of the genotype and allelic distribution of rs2248374 and rs2910686 SNPs among the individuals with AS and the healthy control group, with statistically significant differences (p values = 0.023 and < 0.001, respectively).
Wiśniewski et al., observed no effect of rs2248374 in ERAP2 gene alone on the risk of AS in Polish population.They found associations with ERAP1-ERAP2 haplotypes.The strongest association with AS was observed for rs30187 of ERAP1 gene.The minor T allele and homozygous TT genotype of this SNP significantly increased disease risk (OR = 1.56, 95%CI = 1.22-1.99,p = 0.0004 and OR = 2.52, 95%CI = 1.50-4.25,p = 0.001, respectively) [7].
Paladini et al., also observed no association of ERAP2 rs2248374 with the risk of AS in Sardinian population [26].They found associations with SNP rs75862629. 1 3 vation Funding Authority (STDF) in cooperation with The Egyptian Knowledge Bank (EKB).
are in consistent with the findings reported by Cortes et al. [25], who stated that upon doing a controlled analysis of the ERAP1 interaction with the ERAP2 polymorphism, a significant connection was identified between ERAP2 SNPs and AS instances.Specifically, the rs2910686 SNP exhibited an OR of 1.19 and a p-value of 0.213 × 10 − 2. Also, the study conducted by Ebrazeh et al. [30], revealed a statistical connection between the rs2910686 SNPs and the risk of AS in HLAB27 positive Iranian population.

Limitations of the study
The present research is subject to some limitations, with one notable constraint being the relatively small sample size, which has the potential to impact the validity and generalizability of our findings.Therefore, it is essential to conduct more research with more extensive sample sizes to demonstrate the validity and reliability of our findings.The anticipated focus of future study should be on the therapeutic impacts of these sensitive biomarkers, rather than their diagnostic implications.

Conclusion
The ERAP2 gene SNPs exhibited a significant correlation with AS susceptibility and clinical characteristics.However, it was shown that these SNPs exhibited significant associations with the serum levels of the inflammatory markers including ESR and CRP.Additional investigations of this gene, especially in haplotype analysis, are required to provide a more comprehensive understanding of the genuine contribution of ERAP2 SNPs to the development and progression of AS.

χ 2 :
Chi square test MC: Monte Carlo p: p value for Relation between rs2248374 and different parameters *: Statistically significant at p ≤ 0.05 1 3

Table 1
Comparison between the three studied groups according to demographic data and clinical data HW p1: p value for Chi square for goodness of fit for Hardy-Weinberg equilibrium *: Statistically significant at p ≤ 0.05

Table 2
Univariate regression analysis for alleles of rs2248374 and rs2910686 between cases and controls Group I + II (n = 100) value for Chi square for goodness of fit for Hardy-Weinberg equilibrium p 2 : p value for Univariate regression analysis for comparing with the reference genotype HW p1: p

Table 3
Comparison between the two studied groups according to different parameters *: Statistically significant at p ≤ 0.05

Table 4
Relation between rs2248374 and different parameters in group I

Table 5
Relation between rs2910686 and different parameters in group I

Table 6
Relation between rs2248374 and different parameters in group II

Table 7
Relation between rs2910686 and different parameters in group II